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Abstract. The angular distribution of the emitted electrons, following the two- 
photon ionization of the hydrogen-like ions, is studied within the framework of second 
order perturbation theory and the Dirac equation. Using a density matrix approach, we 
have investigated the effects which arise from the polarization of the incoming light as 
well as from the higher multipoles in the expansion of the electron-photon interaction. 
For medium- and high-Z ions, in particular, the non-dipole contributions give rise to 
a significant change in the angular distribution of the emitted electrons, if compared 
with the electric-dipole approximation. This includes a strong forward emission while, 
in dipole approximation, the electron emission always occurs symmetric with respect 
to the plane which is perpendicular to the photon beam. Detailed computations for 
the dependence of the photoelectron angular distributions on the polarization of the 
incident light are carried out for the ionization of the H, Xe 53+ and U 91+ (hydrogen- 
like) ions. 

PACS numbers: 32.30.Rj, 32.80.Fb, 32.80.Rm 
1. Introduction 

During the last decades, the multi-photon ionization of atoms and ions has been 
widely studied, both experimentally and theoretically. While, however, the majority 
of experiments were first of all concerned with the multi-photon ionization of complex 
atoms, most theoretical investigations instead dealt with the ionization (and excitation) 
of the much simpler hydrogen-like and helium-like systems. For atomic hydrogen, in 
contrast, multi-photon experiments were rather scarce so far (Wolff et al 1988, Rottke 
et al 1990, Antoine et al 1996) mainly because of the lack of sufficiently intensive (and 
coherent) light sources in the UV and EUV region. With the recent progress in the set- 
up of intensive light sources in the EUV and x-ray domain, such as the fourth-generation 
synchrotron facilities or variously proposed free-electron lasers, two- and multi-photon 
studies on the ionization of inner-shell electrons are now becoming more likely to be 
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carried out in the future (Kornberg et al 2002), including case studies on medium-Z 
and high-Z hydrogen-like ions . With increasing charge (and intensity of the light), of 
course, relativistic effects will become important and have been investigated in the past 
for the two-photon excitation and decay (Goldman and Drake 1981, Szymanowski et al 
1997, Santos et al 2001) as well as ionization (Koval et al 2003) of hydrogen-like ions. 
So far, however, all of these studies were focused on the total (excitation or decay) rates 
and ionization cross sections while, to the best of our knowledge, no attempts have been 
made to analyze the effects of relativity on angular resolved studies. 

In this contribution, we explore the angular distribution of the electrons following 
the two-photon ionization of hydrogen-like ions. Second-order perturbation theory, 
based on Dirac's equation, is applied to calculate the two-photon amplitudes including 
the full (relativistic) electron-photon interaction. The angular distribution of the photo- 
electrons are then derived by means of the density matrix theory which has been found 
appropriate for most collision and ionization processes and, in particular, for the - 
nonrelativistic - two-photon angular-dependent studies (Laplanche et al 1986). Since, 
however, the basic concepts of the density matrix theory has been presented elsewhere 
at various places (Blum 1981, Balashov et al 2000), we will restrict ourselves to rather 
a short account of this theory in Subsection 12.11 Apart from a few basic relations, here 
we only show how the angular distribution of the electrons can be traced back to the 
two-photon transition amplitudes. The evaluation of these amplitudes in second-order 
perturbation theory and by means of Coulomb-Green's functions are discussed later 
in Subsections 12.21 and 12 .3| and including the full decomposition of the photon field in 
terms of its multipole components in Subsection 12.41 Using such a decomposition, we 
have calculated the electron angular distributions for the two-photon ionization of the 
Is ground state of hydrogen (H) as well as hydrogen-like xenon (Xe 53+ ) and uranium 
(U 91+ ). By comparing the angular distributions for different nuclear charges Z, we 
were able to analyze both, the effects of the polarization of the - incoming - light 
and the contributions from higher (i.e. non-dipole) multipoles in the decomposition of 
the electron-photon interaction. These results are displayed in Section El and clearly 
show that, with increasing charge Z, the higher multipole components lead to a strong 
emission in forward direction (i.e. parallel to the propagation of the light), while the 
electric-dipole approximation alone gives rise to a symmetric electron emission around 
the polar angle 9 = 90°, similar as obtained by nonrelativistic computations (Zernik 
1964, Lambropoulus 1972, Arnous et al 1973). Finally, a brief summary on the two- 
photon ionization of medium and high-Z ions is given in Section 0] 

2. Theory 

2.1. Density matrix approach 

Within the density matrix theory, the state of a physical system is described in terms of 
so-called statistical (or density) operators (Fano 1957). These operators can be consi- 
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dered to represent, for instance, an ensemble of systems which are — altogether - 
in either a pure quantum state or in a mixture of different states with any degree of 
coherence. Then, the basic idea of the density matrix formalism is to accompany such 
an ensemble through the collision process, starting from a well defined 'initial' state and 
by passing through one or, possibly, several intermediate states until the 'final' state of 
the collision process is attained. 

In the two-photon ionization of hydrogen-like ions, the 'initial' state of the (com- 
bined) system 'ion plus photons' is given by the bound electron \n b j b ^ b ) and the two 
incoming photons, if we assume a zero nuclear spin 1 = 0. For the sake of simplicity, 
we also restrict our treatment to the case that both photons will have equal momentum: 
= k 2 = k, while the spin states of the photons may still differ from each other and 
are characterized in terms of the helicity parameters Ai, A2 = ± 1 (i.e. by means of 
their spin projections onto the direction of propagation k). Of course, the case of equal 
photon momenta k correspond to the most frequent experimental set-up of the two- 
photon ionization of atoms and ions using, for instance, lasers or synchrotron radiation 
sources. With these assumptions in mind, the initial spin state of the overall system 
is determined by the direct product of the statistical operators of the ion and the two 
incident photons 

Pi = Pb ® P 7 <g> p 7 (1) 

or, explicitly, in a representation of the density matrix in terms of the individual 
momenta by 

(n b j b /j b , kAi, kA 2 \pi\ n b j b n' b , k\[, k\' 2 ) 

= (n b j b p b \p b \n b j b fi' b ) (kAi \p 7 \ kAi) (kA 2 |p 7 | kA' 2 ) . (2) 

In the 'final' state of the ionization, after the electron has left the nucleus, we just 
have a free electron with asymptotic momentum p and spin projection m s (as well as 
the bare residual ion with nuclear charge Z). Therefore, the final spin state is described 
by the statistical operator of the emitted (free) electron p e which, in the framework of 
the density matrix theory, can be obtained from the initial-state density operator pi 
owing to the relation (Blum 1981, Balashov et al 2000) 

Pf = Pe = RpiR + ■ (3) 

In this simple relation, R is called the transition operator and must describe the 
interaction of the electron with the (two photons of the) radiation field. Of course, 
the particular form of the transition operator R depends on the framework in which we 
describe the coupling of the radiation field to the atom. As appropriate for high-Z ions, 
below we will always refer to a relativistic treatment of the electron-photon interaction, 
based on Dirac's equation and the minimal coupling of the radiation field (Berestetskii 
et al 1971). 

Instead of applying Eq. in practice, it is often more convenient to re-write 
the statistical operators in a matrix representation. Using, for example, the initial spin 
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Figure 1. Geometry of the two-photon ionization. The photoelectron is emitted 
along the unit vector p = (6, (f>) where 9 is the (polar) angle between the incident 
photon momenta k (chosen as the z-axis) and the electron momentum p. Moreover, 
the (azimuthal) angle <fi defines the angle of p with respect to the x-z plane which, in 
the case of linearly-polarized light, contains the polarization vector u. 



density matrix (J2J), we easily obtain the density matrix of the (finally) emitted electron 
by 

(pm s | p e | pm'J = ^ (n b j b fi b \ p b \ n b j b fx' b ) (k\ 1 | p 7 | kAi) (kA 2 | p 7 | kA' 2 ) 

x M bp (m s ,fi b , Ai, A 2 ) M bp (m' s , n' b , X[, A' 2 ), (4) 
where used is made of the abbreviation 

M bp (m s ,p fe , Ai, A 2 ) = (pm s R k\ 1 , kA 2 , n b j b /j b ) (5) 



in order to represent the transition amplitudes for the two-photon ionization. The final- 
state density matrix (0J) still contains the complete information about the ionization 
process (i.e. the properties of the bare ion and the electron) and, thus, can be used 
to derive all the observable properties of the photoelectrons. Obviously, however, the 
outcome of some considered experiment will depend on the particular set-up and the 
capability of the detectors for resolving the individual properties of the particles. In the 
density matrix theory, this set-up of the experiment is typically described in terms of a 
(so-called) detector operator P which characterizes the detector system as a whole. In 
fact, this detector operator can be considered to project out all those quantum states of 
the final-state system which leads to a 'count' at the detectors; in the language of the 
density matrix, therefore, the probability for an 'event' at the detector is simply given 
by the trace of the detector operator with the density matrix: W = Tr(P p) . 

To determine, for instance, the angular distribution of the emitted (photo-) elec- 
trons, we may assume a detector operator in a given direction p = (9, (j>) [cf. Figure H] 
which is insensitive to the polarization of the electrons 

P = ^2\pm s ) (pm s \, (6) 
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i.e. a projection operator along p and including a summation over the spin state m s 
of the electrons. From this operator, and by taking the trace over the product (P pf) 
with the final-state density matrix (pEJ) . we obtain immediately the electron angular 
distribution in the form 

W(p) = Tr(Pp f ) = — E ( kAl lArl kA i) (kA 2 |p 7 |kA' 2 ) 

Ai(,"i s A1A1A2A2 

x M bp (m s , p 6 , Ai, A 2 ) M fe * p (m s , /x 6 , A' 1; A' 2 ) (7) 

where, for the sake of simplicity, we have assumed that the hydrogen-like ion is initially 
unpolarized. Apart from this additional assumption, however, Eq. flj|) still represents 
the general form of the electron angular distribution for the process of the two-photon 
ionization of hydrogen-like ions. As seen from this equation, the emission of the 
photoelectron will depend on the spin state of the incident photons, defined by the 
photon density matrices (kA|p 7 |kA'). For any further evaluation of this distribution 
function, therefore, we shall first specify these density matrices or, in other words, the 
polarization of the incoming light. For example, if both photons are unpolarized, the 
(two) photon density matrices simply reduce to a constant 1/2, (kA | p 7 | kA') = 5\y/2 
[cf. Appendix, Eq. IA.2j and leads us to the angular distribution 

W^(p) = - E |M 6p (m s ,p 6 ,A 1 ,A 2 )| 2 . (8) 

^ Jb ' ti b m s AiAa 

For many (modern) light sources such as lasers or synchrotron radiation, it is not 
very practical to consider only unpolarized light from the very beginning. In general, 
instead, the angular distribution of the emitted electrons will depend both on the type 
as well as the degree of the polarization of the incident light. For circularly polarized 
light with degree Pc, for instance, the photon density matrix from Eq. (J7J) becomes 
(kA I p 7 I kA') = (1 + XPc) 5\\'/2 and, hence, give rise to the angular distribution 

W™%v) = 1 E ( X + A i p c) (1 + A 2 P C ) |M 6p (m„^,A 1 ,A 2 )| 2 , (9) 

while, for linearly polarized light along the and with a polarization degree Pl, 

the photon density matrix is (kA | p 7 | kA') = 5\\>/2 + (1 — 8\\i)Pl/2. If we evaluate 
Eq. (0) again with this latter density matrix, we obtain the angular distribution 



W £(P) = ~a7oa — TTT E I (!-^) 2 E \M bp (m 8 ,p Jb ,\ 1 ,\ 2 )> 2 
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(10) 



for the electrons as emitted in the two-photon ionization of hydrogen-like ions with 
linearly polarized light. 
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2.2. Two-photon transition amplitude in second- order perturbation theory 

For any further analysis of the electron angular distributions, following the two- 
photon ionization of a hydrogen-like ion, we need to calculate the transition amplitude 
Mb p (m s , fit,, Ai, A2) as seen from Eqs. (|8j)-([rOJ) . This amplitude describes a bound-free 
transition of the electron under the (simultaneous) absorption of two photons. For 
a moderate intensity of the photon field, of course, this amplitude is most simply 
calculated by means of second-order perturbation theory (Laplanche et al 1976) 

M ip (m. )W) A 1) A 2 ) =— ^ ± _ Eb -L- , (11) 

where the transition operator a U\ e* kr describes the (relativistic) electron-photon inter- 
action, the unit vector u\ the polarization of the photons, and where the summation 
runs over the complete one-particle spectrum. In Equation (fTTj) . we added the factor 
VSn^/aE^ in order to ensure that the squared transition amplitude \M bp (m s , fi b , Ai, A 2 )| 2 
has a proper dimension of cross section. From the energy conservation, moreover, 
it follows immediately that the energies of the initial bound state, Eb, and the final 
continuum state, Ef, are related to each other by Ef = Eb + 2E 7 , owing to the energy 
of the incoming photons, which can be written in Hartree atomic units as E^ = k/a. 
Although known for a long time, the relativistic form of the transition amplitude (jllj) has 
been used only recently in studying multi-photon ionization processes and, in particular, 
in order to calculate the total ionization cross sections along the hydrogen isoelectronic 
sequence (Koval et al 2003). In such a relativistic description of the transition amplitude 
(ITTj) . the initial state Vw&M&( r ) = ( r I nbjbUb) and the final state ip pm3 (r) = (r | pm s ) 
are the (analytically) well-known solutions of the Dirac Hamiltonian for a bound and 
continuum electron, respectively (Berestetskii et al 1971). 

As seen from Eq. (J11J) . the evaluation of the transition amplitude requires a 
summation over the discrete (bound) states as well as an integration over the continuum 
of the Dirac Hamiltonian, (i/j u ,E u ). In fact, such a 'summation' over the complete 
spectrum is difficult to carry out explicitly since, in particular the integration over the 
continuum requires the calculation of free-free transitions. This summation, therefore, is 
sometimes restricted to some small — discrete — basis, assuming that the contribution 
from the continuum is negligible. In practice, however, such a limitation seems justified 
only to estimate the behaviour of the cross sections near the resonances where the ion 
is rather likely excited by the first photon into some - real - intermediate state of the 
ion from which it is later ionized by means of a second photon. In the non-resonant 
region of the photon energies, in contrast, the integration over the continuum may give 
rise to a rather remarkable contribution to the total cross section and, hence, has to 
be carried out. Apart from a direct summation over the continuum states, however, it 
is often more favorable to apply Green's functions, at least if these functions can be 
generated efficiently For hydrogen-like ions, for example, such Green's functions are 
known analytically, both in the nonrelativistic as well as relativistic theory (Swainson 
and Drake 1991). 
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2.3. Green's function approach 

As usual, Green's functions are defined as solutions to some inhomogeneous (differential) 
equation 

E-H) G E {r,r') = 5{r - r') (12) 



where, in our present investigation, H refers the Dirac Hamiltonian and E denotes the 
energy of the atom or ion. For realistic systems, of course, such Green's functions are 
not easy to obtain, even if only approximate solutions are needed. However, a formal 
solution is given by (Morse and Feshbach 1953) 

GE(ry) ^\MEmpi : (13) 

including a summation (integration) over the complete spectrum (of H) as discussed in 
the previous section. In the two-photon transition amplitude (fTT|) . therefore, we may 
simply replace this summation by the corresponding Green's function 



M bp (m s ,fi b , Ai, A 2 ) = — —(ip pms {r) au Al e !kr Gf; i+£7 (r, r') au X2 e lkr ' ip nbKbflb (r')) .(U) 



aE^ 

For hydrogen-like ions, the Coulomb-Green's functions from Eq. (|12j) are known 
analytically today in terms of (various) special functions from mathematical physics 
and, in particular, in terms of the confluent hypergeometric function \Fi(a, b; z). Here, 
we will not display these functions explicitly but refer the reader instead to the literature 
(Swainson and Drake 1991, Koval and Fritzsche 2003). For the further evaluation of 
the transition amplitudes (fPfj) let us note only that, also for the one-particle Dirac 
Hamiltonian, the Coulomb-Green's function can be decomposed into a radial and an 
angular part 

„ 1 ^ / gM(r,r') V Km (r) QL(r') -*g|£ (r,r') Km (f) tf_ Km (r') \ 
G E [r,Y) = — - > ,(15) 
rr ™ V ig^(r,r') n_ Km (r) ^L(r') gj£ (r,r>) 0_ Km (r) nl Km (r') J 

where the Q Km (r) denote standard Dirac spinors and where the radial Green's function 
is given in terms of four components g^'ir, r') with T = L, S referring to the large and 
small components of the associated (relativistic) wave functions. The computation of 
the radial Green's function for hydrogen-like ions has been described and implemented 
previously into the Greens library (Koval and Fritzsche 2003); this code has been 
used also for the computation of all transition amplitudes and (angle-differential) cross 
sections as shown and discussed below. 



2.4- Exact relativistic formulation of the two-photon amplitude 

Eq. f!14j) displays the two-photon transition amplitude in terms of the (relativistic) 
wave and Green's functions of hydrogen-like ions. For the further evaluation of this 
amplitude, we need to decompose both, the photon as well as the free-electron wave 
functions into partial waves in order to make later use of the techniques of Racah's 
algebra. As discussed previously for the capture of electrons by bare, high-Z ions 
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(Surzhykov et al 2002a), we first have to decide about a proper quantization axis ( z— 
axis) for this decomposition, depending — of course — on the particular process under 
consideration. For the photoionization of atoms, the only really preferred direction of 
the overall system is given by the photon momenta = k 2 = k which we adopt as the 
quantization axis below. Then, the multipole expansion of the radiation field reads as 

oo 

u A e* kr = u A e^ = £ i L [L]^ ( a£? + «AA« ) , (16) 

L=l 

where [L] = (2L + 1) and the standard notation A^j^ is used for the electric and 
magnetic multipole fields, respectively. Each of these multipoles can be expressed in 
terms of the spherical Bessel functions jiikr) and the vector spherical harmonics T|f A 
of rank L as (Rose 1957): 



Ait = h-i{kr) yj^±T^. x - j L+1 (kr) yj ^yT^ +1 . (17) 

Using the expressions (jlfijl and (|T7|) for the photon field, we can re-write the two- 
photon transition amplitude ()14|) in terms of its electric-magnetic components 

7 L,L'=1 AA' 

x (ip pms \a j\(kr) T X L ) A G Eb+E ^(r, r) a j A >(kr') T% A ,\ip nbKbltb ) » ( 18 ) 
where the coefficients are defined as 

1 if A = L 



e LA = iX Jim ifA = L-l . (19) 



-*Vra ifA = L + 1 

As seen from the expansion (|18j) . we can distinguish between different multipole compo- 
nents such as ElEl, eImI, e1e2, and others owing to the symmetries of the two 
vector spherical harmonics, i.e. due to the particular combination of the summation 
indices L, L' , A, A' in this expansion. In the second line of (|18|). however, the — electro- 
magnetic — multipole matrix elements still contain the wave function i/j pms (r) of the free 
electron with well-defined asymptotic momentum p. In another expansion, therefore, 
we have to decompose it into partial waves to allow for a further simplification of the 
two-photon transition amplitude (|18jl. Again, also the expansion of the free-electron 
wave will depend on the choice of the quantization axis and requires - - by using a 
quantization along the photon momentum — that we have to carry out a rotation of 
the space part of the electron wavefunction from the ^-direction into the p-direction 
(Eichler and Meyerhof 1995) 

VVn s ( r ) = 47r E e ~ iAKf ( l fVf ~ m * l / 2m s I 3fW) 
K f^f 
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where the summation runs over all partial waves /«/ = ±1,±2..., i.e. over all possible 
values of the Dirac angular momentum quantum number Hf = ±(j/ + 1/2) for 
lj = jj ±1/2. In this notation, the (nonrelativistic angular) momentum represents 
the parity of the partial waves and A Kf is the Coulomb phase shift. Moreover, as seen 
from expression (|20j). the partial waves 

^«m.(r)= . 7 (21) 
separate into a radial and an angular parts, where the two radial functions 



->En\ 



P EK (r), g S EK (r) = Q EK (r) 



are often called the large and small components and the corresponding angular parts 
^n^fif) = \lfjfVf) = E m; m s ( l f m i 1 / 2 m s I Jf Hf) Y if mi (r) Xi/2 m a and Q_ Kf/lJ (r) = 

mi(0 Xi/2 m s are the standard Dirac spin 

angular functions. 

Using the partial-wave decomposition (}2~T]) for the free-electron wave function and 
a similar expansion (j!5)l for the Green's functions, we now can carry out the angular 
integration in the transition amplitude ()18j) analytically 

M b (m s ^ b ,X 1 ,X 2 ) = ^^- £ E E * L+L '^P^V A */ 

7 LAL'A' K, f fi f KmTT' 



[L, Lf/ 2 &a & A ' (l L fN - ™ s l/2m s | j //4/ ) 
Kflj [if \a T L \ \ nl m) (kI m\a T L 2 , A , \ K b l b jib 



(P) (22) 

where, apart from the Clebsch-Gordan coefficient (lf[J>f — fn s l/2m s | and some 

constant factors, the angular part of the amplitude is given in terms of the matrix 
elements of the rank L spherical tensor crT^ A = [Y\ ® a]^. These matrix elements can 
be simplified to (Balashov et al 2000) 

K b t[(l b \a T^| KJa'tla^ = [ja, L, A, lf'} 1/2 (j a [l a LM I j b fl b ) 

. f II 1/2 jb ) 

't{0, AO I Zfo) < A 1 L , (23) 

l C' 1/2 Ja J 

by using a proper decomposition in terms of the orbital and spin sub-spaces. The radial 
part of the transition amplitude (fTSj) is contained in (|22|) in the (2-dimensional) integrals 

U^(n f ,K,K b ) = ! gl fKf (r) jA (kr) g T E ^ K (r,r f ) j A ,(kr f ) gjjr') dr dr>, (24) 
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which combines the various (large and small) components of the bound state, the Green's 
function as well as from the free-electron wave. In this notation, again, T = L, S and 
a superscript T refers to the conjugate of T, i.e. T = L for T = S and vice versa. 
In contrast to the angular integrals (J23J), the radial integrals (|2"4")l have to be computed 
numerically. In the present work, all the required integrals for the two-photon transition 
amplitudes (|2*2*Jl are calculated by using the Greens (Koval and Fritzsche, 2003) and 
Racah (Fritzsche et al 2001) programs. 

2.5. Electric dipole approximation 

The transition amplitude ()22j) still describes the full interaction between the electron 
and photon fields. With the explicit summation over all the multipoles of the photon 
field (|16|h it includes the so-called retardation effects or non-dipole contributions. In 
practice, however, the contributions from the higher multipoles decreases very rapidly 
with L and may therefore be neglected; in fact, the computation of these contributions 
also become rather tedious because of difficulties with a stable procedure for the 2- 
dimensional radial integrals (|24j) . In many cases, therefore, it seems justified to restrict 
the summation in (|2*2*jl to just the (dominant) electric dipole term with L = 1 and 
A = L± 1. This 'dipole approximation' is valid if the photon wave length is much larger 
than the size of the atom, i.e. kao <C 1 where cto is the Bohr radius. For the two-photon 
ionization, this condition is well satisfied for most light ions with, say, Z < 30 and for 
photon energies below of the one-photon ionization threshold. 

From the general form (}2"2"|) of the ionization amplitude, the electric dipole appro- 
ximation is obtained by taking L = V = 1 and A, A' = 0, 2 which — owing to the dipole 
selection rules — then also restricts the summation over Kf, i.e. the allowed partial 
waves for the free electron. For the .fT-shell ionization with (completely) circularly 
polarized light, for instance, the final-state electron can only escape in the d%/2 or d§/2 
states. And, as seen from Eq. (j22j) . the dipole transition amplitude is then indeed 
defined by the (second-rank) spherical harmonic, Mb p (m s , A, A) oc Y 2t w - m »+2A(p) 
which (together with Eq.(J£5J)) leads us to the well-known angular distribution 



of the photoelectrons (Lambropoulos 1972, Arnous et al 1973). As expected from the 
axial symmetry of the overall system 'ion plus photons', the angular distribution (|2*5j) 
only depends on 9 but not on the azimuthal angle <fi. For linearly polarized light, in 
contrast, a reaction plane is naturally defined by the photon momentum k and the pola- 
rization vector u and, hence, the axial symmetry is broken. For a linear polarization 
of the incident light, therefore, the angular distribution will depend on both, the polar 
and azimuthal angle and is given by (Zernik 1964, Lambropoulos 1972) 



where the angle = corresponds to an electron emission within the reaction plane [cf. 
Figure 1]. 



W circ {p) = c 4 sin 4 fl 



(25) 



W Vm (p) = 6 + 6 2 sin 2 9 cos 2 + 64 sin 4 9 cos 4 0, 



(26) 
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Figure 2. Angular distributions of the emitted electrons in the two-photon K— 
shell ionization of hydrogen-like ions by means of circularly, linearly and unpolarized 
light. Results are presented for both, the electric dipole ( — ) and the relativistic ( — ) 
approximations and for a two-photon energy which is 40 % above the (one-photon) 
ionization threshold. 



3. Results and discussion 

For the calculation of total two-photon ionization cross sections, the electric dipole 
approximation was recently found sufficient for most of the hydrogen-like ions, and 
not just in the low-Z domain (Koval et al 2003). Even for high-Z ions, for example, 
the total cross sections from the dipole approximation do not differ more than about 
20 % from those of a full relativistic computation, including the contributions from 
all the higher multipoles. Larger deviations, however, can be expected for the angular 
distribution of the emitted electrons which is known to be sensitive to the retardation 
in the electron-photon interaction (Surzhykov et al 2002b). As known, for instance, 
from the radiative recombination of high-Z ions, which is the time-inverse process for 
the one-photon ionization, a significant change in the angle-differential cross sections 
may arise from the higher multipoles and may lead to quite sizeable deviations when 
compared with the dipole approximation (Eichler and Meyerhof 1995). 

In this contribution, therefore, we have analyzed both, the electric dipole and 
the exact relativistic treatment from Eq. (|22|) in order to explore the relativistic and 
retardation effects on the angular distributions of the electrons. Detailed computations 
have been carried out, in particular, for the .ff -shell ionization of (neutral) hydrogen as 
well as hydrogen-like xenon and uranium ions at an energy of both incoming photons 
of Ej = 1.4 • \Ei s \/2 where the E\ s is the (one-photon) ionization threshold. Moreover, 
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Figure 3. Angular distributions of the electrons emitted in the two-photon K- 
shell ionization of the hydrogen-like uranium U 91+ by means of linear polarized light. 
Distributions are shown for the angles = 0°, 45° and 90° with respect to the reaction 
plane; cf. Figure 1. 



to explore the dependence of the relativistic effects on the polarization of the incoming 
light, three cases of the polarization are considered: (i) completely circular polarized, 
(ii) completely linear polarized as well as the case of (iii) unpolarized light. For these 
three ions and types of polarization, Figure 2 displays the angular distributions of 

the electrons as obtained within the dipole approximation ( ) as well as the exact 

relativistic treatment ( — ) which is given by Eqs. (j22j) - (}2"4*)) and includes, therefore, all the 
multipoles in the electron-photon interaction. While, for hydrogen, both approximation 
virtually yields identical results, they start to differ as the nuclear charge Z is increased. 
Instead of a symmetrical emission with respect to the polar angle 9 = 90°, then the 
emission occurs predominantly into forward direction, an effect which is best seen for 
hydrogen-like U 91+ ions. We therefore find, that the non-dipole terms give first of all 
rise to an asymmetrical shift in the angular distribution of the electrons which could be 
observed in experiment. The maxima in the (angle-differential) cross sections, on the 
other hand, are less affected and deviate, even for hydrogen-like uranium, less than a 
factor of 2. 

In Figure 2, all angular distributions are shown as function of the polar angle 9, 
i.e. with respect to the incoming photon beam. As discussed above, this dependence 
of the differential cross sections, |£ = , can be the only one for circular and 

unpolarized light for which the electron emission must be axially symmetric. For linear 
polarized light, in contrast, the emission of the electrons will depend on both, the polar 
angle 9 and the azimuthal angle <fi. For this polarization, Figure 2 only displays the 
angular distributions within the reaction plane, i.e. at <fi = 0°. To explore, in addition, 
also the ^-dependence of the two-photon ionization by linear polarized light explicitly, 
Figure 3 shows the corresponding angular distributions ^(9, (p) for the three particular 
angles <fi = 0°, 45° and 90° with respect to the reaction plane; here, the left inlet (<fi 
= 0°) is the same as shown in Figure 2 in the middle column for U 91+ ions. Again, 
the results from the electric dipole approximation are compared with those from a fully 
relativistic computation. As seen from Figure 3, the most pronounced effect of the 
higher multipoles arise for an electron emission in a plane, which is perpendicular to the 
photon polarization vector (<f) = 90°). In such a — perpendicular — geometry of the 
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Figure 4. Angular distributions of the electrons emitted in the two-photon if -shell 
ionization of the hydrogen-like uranium U 91+ by circular polarized light with different 
degrees of polarizations Pc = 0, 0.5, 0.7. and 1. 

experiment, the cross sections from the exact treatment show strong forward emission 
of the photoelectrons while the electric dipole approximation (|26|) . in contrast, results 
in a completely isotropic emission, if seen as function of the polar angle 9. 

Until now, we considered the two-photon ionization of hydrogen-like ions by 
either completely polarized (linear: Pl = 1; circular Pc = 1) or unpolarized light 
(Pl = Pq = 0). In most experimental investigations on two- (and multi-) photon 
processes, however, the incident radiation is typically polarized with some given degree 
of polarization < Pc, -Pl < 1. Apart from the type of the polarization of the in- 
coming light, therefore, we shall study also how the angular distributions depend on 
the degree of polarization. Figure 4, for instance, displays the angular distribution from 
the K— shell ionization of hydrogen-like U 91+ ions by means of circular polarized light 
with a degree of polarization Pc = 0.0 (unpolarized case), 0.5, 0.9 and 1.0. As seen 
from this figure, the probability for an electron emission increases at angles around 9 
= 60° but decreases (towards zero) in forward and backward direction as the degree 
of polarization is increased. In particular the behaviour near 9 = 0° and 9 = 180° 
can be easily explained if we consider the conservation of momentum in the overall 
system. Since, for completely circularly polarized light, the (total) spin projection of 
photons on the quantization axis (which is chosen along the photon momenta k) becomes 
Ai + A2 = ± 2, it obviously can not be compensated - in the final state - if the electron 
is emitted parallel (or antiparallel) to the incoming light and hence its spin projection 
is jJLf — m s = ± 1/2. For unpolarized light, in contrast, the photons may have different 
helicities and, therefore, the projection of their angular momentum Ai + A2 = may be 
conserved under a forward and backward non spin-flip electron emission. 
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4. Summary 

In this paper, the two-photon ionization of hydrogen-like ions has been studied in the 
framework of second-order perturbation theory and the relativistic description of the 
electron and photon fields. That is, exact Dirac bound and continuum wave functions 
were applied for the description of the electron to reveal the importance of relativity 
on the angular distributions of the emitted electrons. Moreover, relativistic Coulomb- 
Green's functions are used to perform the summation over the complete Dirac spectrum 
as needed in second-order perturbation theory. 

To understand the angular distributions of the emitted photoelectron and, in 
particular, the influence of the polarization of the light on this emission, density matrix 
theory has been utilized to 'combine' the two-photon transition amplitudes in a proper 
way. Calculations are carried out for the i^-shell ionization of the three (hydrogen-like) 
ions H, Xe 53+ and U 91+ . From the angular distribution of the electrons for different 
types (linear, circular, unpolarized) and degrees of polarization (i.e. in going from the 
completely polarized to unpolarized light), it is clearly seen that the angular emission 
depends much more sensitive on the contributions from higher multipoles than the total 
cross sections. Two rather pronounced effects, for example, concern the (asymmetrical) 
forward emission of the electrons as well as a significant change in the electron emission 
for linear polarized light, if the electrons are observed perpendicular to the reaction plane 
[cf. Figure 4]. Both effects are enhanced if the nuclear charge of the ions is increased. 

An even stronger influence from the non-dipole terms (of the radiation field) is 
expected for the spin-polarization of the photoelectrons. Similar as in the present 
investigation, density matrix theory provides a very suitable tool for such polarization 
studies. A detailed analysis of the polarization of the photoelectrons, emitted in the 
two-photon ionization of hydrogen-like ions, is currently under work. 
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Appendix A. Photon spin density matrix 

A pure (i.e. completely polarized) state of the photon can be characterized in terms 
of a polarization unit vector u which always points perpendicular to the (asymptotic) 
photon momentum k. Of course, this polarization vector, u, can be re-written by means 
of any two (linear independent) basis vectors such as the circular polarization vectors 
u±i which are (also) perpendicular to the wave vector k and which, for u +i respective 
u_!, are associated with right- and left-circular polarized photons (Blum 81). In such 
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a basis, the unit vector for the linear polarization of the light can be written as 

u(x) = ( e~* u +1 + e* u_ x ) , (A.1) 

where x is the angle between u(x) and the x— 2 plane. 

While a description of the polarization of the light in terms of either the circular 
polarization vectors u±i or the linear polarization vector (jA.lj) is appropriate for 
completely polarized light, it is not sufficient to deal with an ensemble of photons which 
have different polarization. Such a — mixed — state of the light is then better described 
in terms of the spin-density matrix. Since the photon (with spin 5 = 1) has only two 
allowed spin (or helicity) states |kA) , A = ±1, the spin-density matrix of the photon 
is a 2 x 2 matrix and, hence, can be parameterized by three (real) parameters: 



1 



1 + P c Pl e~ 2ix 

kA |Ar|kA') = i [ |, (A.2) 

2 I P L e 2i * 1 - P c 



where < Pl < 1 and —1 < Pc < 1 denote the degrees of linear and circular 
polarization, respectively. The angle Xi moreover, represents the direction of the 
maximal linear polarization of light. 

Of course, the choice of the parameters Pl, Pc and \ is not unique and many other 
- equivalent — sets of three real parameters could be applied to characterize the photon 
spin density matrix (|A.2|I . In the analysis of experimental data, for instance, one often 
uses the three Stokes parameters to describe the polarization of radiation. The Stokes 
parameters can easily be expressed in terms of the (two) degrees of polarization, Pl and 
Pc, and the angle x as: 

P 1 = P L cos 2 X , P2 = Pl sin 2 X , P 3 = Pc • (A.3) 

The use of the Stokes parameters leads to the familiar form of the spin density matrix 
(Blum 1981, Balashov et al 2000) 

P 3 Pi - iP 2 

kA I p 7 1 kA' ) = ^ I _ i ^ I (A.4) 
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